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1 Introduction 



The dynamics of particles in curved spacetime has been an important part of the theories of 
gravitation. The dynamics of test particles, i.e., particles without any internal structure was 
studied in the early stages of the development of general relativity. Recently, there has been a 
growing interest in the problem of motion of spinning particles (some references may be found 
in [§). 

Rietdijk and van Holten [||] has obtained a general set of equations of the symmetries of 
the spinning particles in an arbitrary curved spacetime. Similar equations has been obtained 
in from a more general action. These equations constitutes a Grassmann valued extension 
of the Killing equations for the symmetries of the spacetime manifold. These ideas was fur- 
ther investigated in [Q] and was applied to the case of Schwarzschild spacetime in [Q] and to 
Reissner-Nordstrom spacetime in The canonical structure of the theory and the equivalent 
Dirac formulation of the quantum theory has been derived in Equations for the world-line 
deviations of spinning particles was obtained in [Js[] . 

The aim of the present paper is to study the motion of spinning particles in the spacetime of an 
exact plane gravitational wave following the ideas presented in the references mentioned above. 
This problem has been studied previously in the framework of Papapetrou-Dixon description 
of spinning particles in [|l|] and [g, |l0[ |ll|]. In [[I]] a linearized version of Papapetrou-Dixon 
equations was applied to the case of a weak plane gravitational wave of arbitrary polarization 



and harmonic profile. In QlOp a linearly polarized wave of arbitrary profile (and in particular, a 
square gravitational wave pulse) was considered and a set of exact solutions was presented in 
some rather special circumstances. The general solution of those equations could only be found 
numerically, of course. 

There is a simpler version of Papapetrou-Dixon equations. This is basically obtained by 
replacing the usual supplementary condition p • s = by the so called Pirani supplementary 
condition v • s = 0. A similar calculation has been performed along the lines of [|lj|, but using 



Pirani condition instead [|12|]. It led to conclusions very close to those of [jlj]. 

In this paper we find the general solution for motion of spinning particles in the spacetime of 
an exact plane gravitational wave. We use only two of the generic constants of motion, namely 
the Hamiltonian and the supercharge. The results may be considered as a generalized alternative 
to the results of [[Ij] and [10]. The solutions we present here, give more physical insight into the 



behaviour of spinning particles and may be of interest in the important problem of the detection 
of gravitational waves. In the subsequent sections we first give brief reviews of spinning particles 
and gravitational waves and then solve the equations of motion exactly, next we study the world- 
line deviations of spinning particles and finally present our physical conclusions. An appendix is 
devoted to the small-/i limit of the particle's trajectory and the evolution of its spin. 

2 Spinning particles 

The equations of motion of a classical fermion may be obtained from the action []7| ^ 

S = m f dT(\g^(x)x»x v + J- g ( x )^^-), (1) 
J 2 2m ut 



x^ being commuting and ip^ anticommuting variables with /i, v, ■ ■ ■ running over the dimensions 
of the spacetime. The above action will lead to the following equations of motion 



n2 u • 

-j^ = —Wvxpx^r, (2) 



Dt 



0, (3) 



where ^ := + T» Xv x x x u , BgL ■- ^ + x x T» Xu ^, a dot denotes and R\ Xfl := 
d x T^ - c^ A + T^r^ - T^ p T p KX is the Riemann curvature tensor. 

There are four types of generic (i.e. existing in any spacetime) conserved quantities connected 
with the action ([!]). These are the Hamiltonian H ,the supercharge Q, the dual supercharge Q*, 
and the chiral charge T+, given respectively by 

I 777 

H = —g^(x)p^ u = -g^(x)± fl ± u , (4) 

Q = p^V, (5) 

Q* = ^V^ge^xP^^, (6) 



3! 

r* = j^^u^Vr^, (?) 

The relativistic spin of the particle is described by the antisymmetric tensor 

sT = (8) 

whose spacelike components represent the particle's magnetic dipole moment and the timelike 
components s° l correspond to the electric dipole moment. For free Dirac particles we require the 
dipole moment to vanish in the rest frame. Hence we set 

9flX {x)x x s^ = 0, (9) 

or equivalently 

g^(x)x^ u = 0, (10) 

That is Q = 0. The conservation of Q guarantees that this condition can be satisfied at all times 
irrespective of the presence of external fields. 

We may recast Eqs. (|2|) and (|3j> in terms of s^ v as follows 

Dt 2 - 2m K vx p x a » (ID 
Ds^ 

~Dt~ = °' <12) 

where := s pu + T^ x x K s Xu + T v kX x k s^ x . The particle's spin may also be represented via a 
spin four-vector given by 

:= -^=e^ KX x u s KX . (13) 
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In general relativity, only relative acceleration has an observer- independent meaning. The rel- 
ative acceleration of two nearby particles can be determined by measuring the rate at which the 
particles world-lines deviate from each other. The world-line deviation equations are given by 



V 



R^ KXv x K x v n x s ab R^ 



uab~ 



Dt 2 2m Dt 



Dt 



-L(s ab V x R» uab i»n x + J ah R» vab x v ), (14) 
2m 

s K ^R u K a pn a x?, (15) 



where V x R^ UKp = d x R% Kp + T^ xp R^ VKp - T? vX R^ Kp - T^ KX R^ p - pX R^ K p and 
n M , J ab are defined for a one-parameter congruence of solutions (x m (t; X),ip(r; A)) to (|J) and 
(^) as follows 

dx^ b Ds ab 

U ~ dX ' ~~ DX 

Equations ([T4])-([T5|) are generalization of the well-known equation of geodesic deviation. The 
first term in the r.h.s of (Q represents the usual effect of the spacetime curvature on the devia- 
tions of world-lines, and the other terms are due to the spin-curvature coupling. These equations 
may be derived by a procedure similar to the derivation of the equation of geodesic deviation, or 
by alternative methods [§]. 

In the next section we seek solutions to Eqs. (|), (0), (|I|) and®,© in a gravitational 
wave spacetime. 



3 Plane gravitational waves 



The most general plane-fronted parallel-rays gravitational wave is given by the metric [13] 

ds 2 = -dudv + F(u, x a )du 2 + dx a dx a , (16) 

where (u,v) are null coordinates, x a are spacelike transverse coordinates, and F(u,x a ) is a 
solution of the two-dimensional Laplace equation with respect to x a . Plane gravitational waves 
are the particular case where F is quadratic in x a . Here we consider a wave given by the metric 

ds 2 = —dudv — K(u, x, y)du 2 + dx 2 + dy 2 , (17) 

where K(u,x,y) = f + (u)(x 2 — y 2 ) + 2f x (u)xy. The arbitrary functions / + and f x corre- 
sponding to two independent polarisations of the wave. We consider a wave of harmonic profile 
(suitably chosen to reduce to the metric given in \^ in the weak field limit, in the so called group 
coordinates) 

f+(u) = huj 2 sin(um), (18) 
/ x (u) = huJ 2 cos(wu), (19) 

h being the dimensionless wave amplitude. 
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4 The particle's trajectory and spin 

1F71) and suitable initial 



In this section we seek solutions to Eqs. (Q) and (|TT])- ( JT3| ) together with (] 

,v,x,y respectively. From Eq. ( 1 

0, (20) 



conditions, /x = 1, 2, 3, 4 correspond to u, v, x, y respectively. From Eq. (12|) we have 

ds 13 



dr 
ds u 



0, (21) 



dr 

ds u 1 . , A 1 



cir 2 2 

The three remaining may be found from Eq. d9|) with v = 1, 3, 4 



-i^tts 14 + -i-^us 13 . (22) 



us 12 = 2is 13 + 2ys 14 , (23) 
us 23 = -(i; + 2i^n)s 13 -2ys 34 , (24) 
us 24 = -(i, + 2Kii)s u + 2xs u ,. (25) 



From Eq. (Mil) we obtain 



cir 2 

d 2 x 1 . 13 1 

dr 2 " 2m ^ 2 

cir 2 " 2m raUS 2 



(26) 



— i^us 13 - -l^n 2 , (27) 



— i^ ra ns 14 - -i^xi 2 , . (28) 

m 



The other component of the trajectory may be found by considering H = — y from which, Eq. 
(Q) results in 

u = -(1- K + x 2 + y 2 ), (29) 



We fix the gauge such that the solution to Eq. (26) takes the following form 



u = t. (30) 

Also Eqs. (|20"|)-(|2"l|) result respectively in 

s 13 (r) = constant = a, (31) 

s (t) = constant = (3. (32) 



So Eqs. (|27|)-(|28D become 



x = — hoo 2 sin(cdT) — huj 2 s\u.(ujt)x — hoj 2 cos(ur)y, (33) 
m 

R 

y = hu 2 sinCwr) — hu; 2 cos(urr)x + hoo 2 cos(ojt)x/, (34) 

m 

respectively. With £ = x + iy and c = a — i(3 the solution to the above system of coupled 
differential equations is 

U J ~ 2m m(4-/i 2 ) + 2m(4 - h 2 ) ^ P 

+ 7T^^- l ( 1+a )^ + ^ + -^ C 7e-( 1+fe )^, (35) 
{l + a) z {l + by 



where a 



4h w - t | 1 ^ q constants determined by initial conditions, and 



overbar denoting complex conjugation. For h < I both x and y components contain sin and 
cos terms only, resulting in a combined oscillation in the x — y plane. For h > | exponential 
terms appear in both x and y. In this case oscillatory terms are modulated by either exponentially 
decreasing or increasing terms and the resulting motion differs significantly from the previous 
type. We do not investigate this type of motion here. 

It is also interesting to find expressions for the Grassmannian variables rp^. Using Eqs. @, 



(10), and (17) we obtain 



1>\t) 



a constant Grassmannian 

^g(r) 



- + — fir) + — gij)) 
u u u 



(36) 

(37) 

(38) 
(39) 



where f(r) and g(r) are solutions to df(r)/dT = —\j1K x u and dg{r)jdT = —l/2K y u re- 
spectively. The immediate result of these expressions is the vanishing of the dual supercharge 
and the chiral charge: 

Q* = r\ = o 

Having the exact solution to Eqs. (|20|)-([2l|) and ([£(])-( 28), it is now straightforward to find solu- 
tions to Eqs. (^2|)-( 25 ) and (|29|). So we have found in principal the exact solutions to the equa- 
tions for the trajectory and the spin of the particle. The resulting expressions are algebraically 
involved. To get physical insight into the nature of the resulting solutions, one may use a weak 
field limit of the solutions. We have gathered in the appendix some of the linearized solutions. 



5 The world-line deviations 



From (17) and (15 ) we obtain 



dJ 13 _ dJ lA 
dr dr 



whose solutions are 



J 13 (r) 
J 14 fr) 



constant = e, 
constant = 5. 



Equations (|14j) results in 



'2„1 



d n 

One may choose the following solution to the above equation 

n x (r) = 

Now inserting the above relations in equations ( Jl~4| ) we get 



d 2 n 6 {r) 

dr 2 
d 2 n 4 {T) 

dr 2 



+hu> 2 sin(<jjr)n 3 (r) H huj 2 sin(a;r) 

m 

—hoj 2 sin (wr)n 4 (r) huj 2 sin(a;r) 

m 



(40) 
(41) 

(42) 
(43) 

(44) 
(45) 
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The other component, n 2 (r) may be obtained from the condition x ■ T j L = 0. An approximate 
solution to (Q and @) is 

n 3 (r) = A + h((— - A) (lot -sin( lot))), (46) 
m 

n 4 (r) = B + h(( + B)(wr-sin(wT))), (47) 

m 

in which A, B denote the initial separations of the two spinning particles at rest near the origin of 
the coordinates with the same -j- and These show how the relative accelerations of spinning 
particles depend on the spin. The terms independent of the spin are due to the usual tidal force. 



6 Discussion 



The projection of the particle's trajectory onto the x — y plane is given by (35). This is depicted 
for a particle initially at rest at the origin in Fig.|] for some set of the parameters. As the figure 
shows the path is an ellipse with a slow precession. The dimensions of this ellipse depend on the 
particles's spin. For a given frequency, this precession gets slower for smaller values of h. 

Using a weak field approximation (see the appendix) give some insight into the nature of the 
exact solution we found. Let us examine the trajectory and spin of the particle in some rather 
special situations. For a particle at rest at the origin of coordinates at r = (which we consider as 
the time in which the wave arrives), if the initial orientation of the spin is chosen to be a = (3 = 0, 
and 7 / then 

(»(r) l «(T) I i(r) ) ii(r)) = (r, r, 0, 0), (48) 
that is the particle remains at rest at the initial coordinates, and 



(s 1 (t),s 2 (t),s 3 (t),s 4 (t)) = (-7,7,0,0), 



(49) 



its spin remains unchanged. These are in agreement with the results of [jl|] and [10] 
For 7 = and a ^ 0, / 0, the particle's trajectory is 



x(t) 

y( T ) 

v(t) 



ltd!. , , OL , , Gi , . 

sm wr H smlhuT) cos(/mjt), 

2m m m 2m 

P , h/3 . P . P 
1 sin(wr) sm(rajT) — - — cos(no'r), 



2m m 

T. 



in 



2m 



(50) 

(51) 
(52) 



and its spin is 



s 2 (r) 



s 3 (r) 
At) 



2huj 

m 
2hu 

m 

P, 

—a. 



aP(l 
aP(l 



cos(wt)), 
cos(wt)), 



(53) 

(54) 

(55) 
(56) 



It is now clear from Eqs. (5C)-(52) that the particle's displacement along the direction of the wave 
propagation (z = ■ ! ^) is much smaller (0(h 2 )) than its displacement in the transverse plane, 
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even though this transverse displacement is practically very small (0(h)) itself. This is similar 
to the behaviour of spinless particles. 

For general initial spin orientation (a / 0,/3 / 0,7 / 0), the transverse components of 
the particle's spin have oscillations of very small amplitudes (of the order of h) around their 
initial orientation. The component of the spin parallel to the direction of the wave propagation 
s z = s T, s remains unaffected, again in agreement with [p. 



The relative acceleration of two nearby spinning particles is given by @ and ©. They 
show that the relative acceleration oscillates with the frequency of the wave. 

As argued in [14] we may recover Eqs. @ and ( [TT| ) -(|l2|) from the so called Papapetrou- 
Dixon equations by neglecting the terms quadratic in spin. On the other hand, as we saw in 
our calculations above, the contribution from spin are usually small. Thus it seems neglecting 
the very small contributions from terms quadratic in spin makes no significant loss of physical 
content. Hence in the lack of exact solutions to Papapetrou-Dixon equations in gravitational wave 
spacetime, our results may be generalized to them as well (at least in spacetime regions far from 
the strong sources of these waves). 
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Figure 1: 

The particle's trajectory projected onto transverse plane for — 

< r < 100000. 
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appendix 

In this appendix we gather some linearized solutions for the particles trajectory and spin. Putting 
C(0) = C(0) = 0, Eq. (35) reduces to first order in h to 



£( r ) — + !^£ e iwT _ l 2^ e -i^r _ c + 2ic 
2m 2m 2m 4m 



e ihuiT _|_ ^ic c c ~ihu)T 

4m 



or equivalently 



ot ho. a. 
x(t) = sin(a;r) H sin(/io;r) — ^— cos(hujT), 



y( T ) 



P ,h/3 p 

1 smiur) sm nwr j 

2m m m 



a 

2m 
_P_ 
2m 



COs(flLdT). 



Thus from Eq. ( |29| ) we obtain (with v(0) = 0) 

v( T ) = r + 0(h 2 ). 



Eq. (30) and ft 



determine the particle's trajectory. FromEq. (22) we obtain 
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(57) 

(58) 
(59) 

(60) 
(61) 



where 7 is a constant of integration. 



We may now put Eqs. (31)-(32) and (plj) into (E3h-(25) to find the other components of s^. 



The results may be fed into (|13|) (with e 



1234 



+1) to give s A 



2huj 

m 
2huj 

m 



a(3(l — cos(wt)) — 7, 

a(3(l — cos(wt)) + 7, 

(3 07(1 — cosfun - )), 

m 

—a — — /?7(1 — cos(wt)), . 
m 



(62) 
(63) 
(64) 
(65) 
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